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ABSTRACT: The Klein-Gordon equation is a useful test arena for quantum cosmolog- 



ical models described by the Wheeler-DeWitt equation. We use the decoherent histories 
approach to quantum theory to obtain the probability that a free relativistic particle crosses 
a section of spacelike surface. The decoherence functional is constructed using path in- 
tegral methods with initial states attached using the (positive definite) "induced" inner 
V_i ! product between solutions to the constraint equation. The construction is complicated by 

bX). 

the fact that the amplitudes (class operators) calculated using a path integral typically do 
not satisfy the constraint equation everywhere, but we show how they may be systemat- 
ically modified in such a way that they do satisfy the constraint. The notion of crossing 
a spacelike surface requires some attention, given that the paths in the path integral may 
cross such a surface many times, but we show that first and last crossings are in essence 
the only useful possibilities. Different possible results for the probabilities are obtained, 
depending on how the relativistic particle is quantized (using the Klein-Gordon equation, 
or its square root, with the associated Newton- Wigner states). In the Klein-Gordon quan- 
tization, the decoherence is only approximate, due to the fact that the paths in the path 
integral may go backwards and forwards in time. We compare with the results obtained 
using operators which commute with the constraint (the "evolving constants" method). 



1. INTRODUCTION 



If the state of a quantum system obeys a wave equation of the form 



= 0, 



(1.1) 



how do we extract probabilities from the wave function? Quantum cosmological models 
are described by precisely such a wave equation - the Wheeler-DeWitt equation - where 
H is the total Hamiltonian of the matter and gravitational fields [1,2]. A simpler but 
important example of an equation of this type is the Klein-Gordon equation 



The Wheeler-DeWitt equation for simple models has the form of a Klein-Gordon equa- 
tion in a general curved spacetime background with a spacetime dependent mass term. 
Traditional approaches to relativistic quantum theory note the various difficulties of inter- 
preting the Klein-Gordon equation and then pass quickly on to quantum field theory. The 
Wheeler-DeWitt equation, however, in its full form, already represents a second-quantized 
field theory. To work with the Wheeler-DeWitt equation we must therefore return to wave 
equations of the Klein-Gordon type and understand how to overcome their difficulties with- 
out resorting to second quantization. These questions are closely related to the problem 
of time in quantum gravity [3,4,5]. Furthermore, it is likely that they will also be present 
in other approaches to quantum gravity, such as the loop variable approach [6]. 

A number of approaches to this problem have been considered. One is the decoher- 
ent histories approach, in which the wave function is associated with a set of histories to 
which probabilities may be assigned [7,8,9]. Another approach involves the use of operators 
which commute with H [10,11,12,13,14,15]. The purpose of this paper is to explore the 
application of these methods to some simple probability questions involving the relativis- 
tic particle in flat spacetime. We will concentrate on the decoherent histories approach, 




(1.2) 
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extending and developing earlier work in this area [1,16]. At the present (rather early) 
stage of development of this field it is not possible to say whether it is equivalent to an 
operator approach, and an important part of this paper is to compare the approaches 
where possible. Indeed, there are many different and potentially inequivalent quantization 
methods applicable to simple parametrized systems [17,18]. 

We will focus on the following question: given a solution to the Klein-Gordon, what 
is the probability of finding the particle in a spatial region A of a spacelike surface? 
The question is clearly a simple one, but it turns out to expose some subtle aspects of 
the decoherent histories approach applied to parameterized systems, and is an important 
test of the formalism in a familiar situation. Furthermore, it may also be regarded as a 
preparatory exercise for the treatment of more complicated quantum cosmological models, 
which will be considered elsewhere. 

We begin by briefly reviewing the various aspects of the formalism relating to the wave 
equations (1.1), (1.2). 



The inner product traditionally associated with the Klein-Gordon and similar equa- 
tions is the Klein-Gordon inner product, 



It is evaluated on a spacelike surface E, and is independent of the choice of such surface if 
\& and $ are solutions to the Klein-Gordon equation. This inner product is, however, not 
positive definite. When a separation into positive and negative frequencies is possible, it 
is positive on the positive frequency sector and negative on the negative frequency sector. 



1(A). Inner Products 




(1.3) 
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An essential part of our approach here is to use the methods of refined algebraic 
quantization, in which we work with the so-called induced (or Rieffel) inner product [16,19]. 
This starts through the introduction of an auxiliary inner product 

(V,$) A = J d 4 x V*(x)$(x) (1.4) 

We then consider eigenstates of the constraint 

H^ Xk = \^ Xk (1.5) 

where k is a degeneracy label. These are normalizable in the auxiliary inner product via 

^Xk^X'k')A = S(X-X')8(k-k') (1.6) 

The induced inner product between solutions to the constraint then consists of dropping 
the 5(X — A') term on the right and taking the limit A, A' — > 0. This produces a well-defined 
positive definite inner product on solutions to the constraint. This procedure may also be 
understood by starting from the observation that solutions to the constraint (1.5) may be 
written in the form S(H — X)\x) for some fiducial state |x). The induced inner product 
between such states is then effectively equivalent to replacing [5(H — A)] 2 with 5{H — A), 
that is, taking 5(H — A) to be a projection operator (which it clearly is in the case of a 
discrete spectrum). 

When a split into positive and negative frequency solutions is possible, ^ = + \f ' , 
the induced inner product coincides with the Klein-Gordon inner product but with the 
sign of the negative frequency sector changed so as to make the product positive: 

0/ d> = (VJT, $)j = (VJT+, § + ) KG - (*-,*-) KG (1.7) 

For the free relativistic particle, we could quite simply have defined an inner product by the 
object on the right, and this is well-defined since the positive and negative frequency sectors 
do not interact in this case. The advantage of the induced inner product, however, is that 
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provides a good inner product even when the split into positive and negative frequencies 
is not possible. (See Refs. [16,19] for more details). 



1(B). An Operator Approach 



Given the constraint equation (1.1) and its inner product structure we would like to 
be able to assign probabilities to various dynamical variables of interest. It is generally be- 
lieved that the interesting dynamical variables are those that commute with the constraint 
H [10,11,12,13,14,15]. This is because the constraint is associated with reparametriza- 
tion invariance (diffeomorphism invariance in the general case) and we are interested in 
variables that are invariant. 

Because the wave equation is not of the Schrodinger type, it does not have an external 
time variable, so we cannot talk about the value of a variable at a particular "time". 
Instead, "time" is somehow encoded in the variables already present in the wave equation. 
In the Klein-Gordon equation, for example, we might be interested in the value of the 
spatial coordinate x at given x°. We might equally be interested in the value of x°, say, 
at given x 1 . We need operators commuting with H which express these quantities. 

Suppose we are interested in the operator corresponding to the value of A when B 
takes the value r. The appropriate operator is 



ordering is chosen in this expression, although note that it is not always possible to make 
it self-adjoint). It is readily verified that this operator commutes with H. The study 
of operators of this type is the basis of the operator approach (sometimes know as the 
"evolving constants" method). The spectrum of the operator Eq(1.8) is computed from 
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(1.8) 



where A(s) = e lHs Ae lHs , and similarly for B(s) [14,15]. (W 



e assume a suitable operator 



which one may compute a projection operator P a say, onto a range of the spectrum. The 
associated probability is then of the form Tr(P a p). 

1(C). Decoherent Histories 

In this paper we are primarily concerned with the decoherent histories approach to 
quantum theory, and this provides a second method of calculating probabilities of interest 
On the face of it this method is quite different to the operator method outlined above, and 
as stated in the Introduction, there is no particular reason to assume that the methods are 
equivalent. 

In the usual formulation in non-relativistic quantum mechanics [7,8,9], the central 
object of interest is the decoherence functional, 

D{a,oi) = ^r[p f C a pCl) (1.9) 
where the histories are characterized by the class operators C a , which satisfy 

J2 C * = 1 ( L1 °) 

a 

In the simplest case, for non-relativistic systems, the class operators are given by time- 
ordered sequences of projection operators 

C a = Pa n (tn) ■ ■ ■ PaAh) (1-H) 

where a denotes the string of alternatives a\, oli • • ■ a n . The theory is, however, more 
general than this and we will exploit this generality here [1,20]. We have included the pos- 
sibility of both an initial state p and a final state p j (normally taken to be proportional to 
the identity), and we therefore have to include the normalization factor N = (Tr(pj-p)) -1 . 

Intuitively, the decoherence functional is a measure of the interference between pairs 
of histories a, a'. When its real part is zero for all pairs of histories with a ^ a/, we say 
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that the histories are consistent and probabilities 



p(a) = D(a, a) 



(1.12) 



obeying the usual probability sum rules may be assigned to them. Typical physical mech- 
anisms which produce this situation usually cause both the real and imaginary part of 
D{ql,ol) to vanish. This condition is usually called decoherence of histories, and is re- 
lated to the existence of so-called generalized records [7,21]. Note also that when there is 
decoherence, using (1.10) the probabilities may be written 



In its application to the parametrized systems considered here, the initial and final 
states are attached to the class operators using the induced inner product scheme [16]. So 
for pure initial and final states, the decoherence functional is 



where o here denotes the induced inner product. (We may of course sum over initial or 
final pure states to get mixed ones.) The key question for the models considered here is 
then the construction of class operators corresponding to questions of interest. This is the 
step analogous to the construction of operators above in Eq.(1.8), and therefore the class 
operators must somehow incorporate reparametrization invariance. It should be stated at 
this stage that there does not at present seem to be a completely clear and unambiguous 
prescription for constructing the class operators, and part of the aim of this paper is 
therefore to explore possible constructions and examine their properties. 

One very natural approach to calculating the class operators for reparametrization- 
invariant systems is to use path integrals [16,1]. For the relativistic particle these have the 
form 




(1.13) 




(1.14) 




(1.15) 
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where 

r-T 



g a (x",T\x',0) = J Vx^expl-iJ 



ds 



+ m 



(1.16) 
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The path integral for g has the form of a non-relativistic propagator. The sum is over paths 
from spacetime points x' to x" where the paths are restricted in some way defined by the 
coarse-graining a. For example, we might be interested in the probability that the particle 
passes through some region of spacetime, or not. More details about the construction 
of this object, including the specification of the range of the T integration, will be given 
below. 

An issue which arises with this definition is that these class operators often do not 
everywhere satisfy the constraint equation with respect to their end-points. As we shall see 
in particular examples, they satisfy the constraint except on the boundaries of the regions 
defining the coarse grainings. This is an issue because the induced inner product is only 
defined between solutions to the constraint. The auxiliary inner product is still defined, 
but in operating with such a class operator on the initial state, we are stepping out of 
the constraint surface. Fortunately, we will find in particular examples that the problem 
is easily fixed by a small amount of intuitively sensible doctoring on the class operators, 
guided by the requirement that path integral methods agree with operator methods. In 
particular, we will follow the suggestion of Hartle and Marolf, which, loosely speaking, is 
to replace C a by a new object C' a which satisfies the constraint equation everywhere, and 
satisfies the same essential boundary conditions as the path integral-defined object C a [16]. 

It is easily seen that the operator and decoherent histories approaches are different, 
with no guarantee of their equivalence in general. They both deal with objects which are 
compatible with the constraint (the operator (1.8) and the class operators C a ), but the 
operator method looks at a projection operator onto ranges of the spectrum of (1.8), whilst 
the decoherent histories approach works with the class operators C a and the decoherence 
functional (1.14). The class operators are not projections in general, and in some sense 
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are generalizations of projections to non-commuting alternatives, so the two formalisms 
are quite different. However, it is known that when there is exact decoherence in the 
decoherent histories approach, the probabilities for histories Eq.(1.12) may be written in 
the form Tr(R a p) where R a is a projection operator (corresponding to the existence records 
mentioned earlier [7,21]). 

1(D). Summary of this Paper 

As stated, the main aim of this paper is to derive expressions for the probability 
of crossing a spacelike surface in relativistic quantum mechanics, using the decoherent 
histories approach, and the similar operator approaches. 

To get a feel for the formalism for reparametrization-invariant systems, we start in 
Section 2 by applying the formalism to the non-relativistic particle in parametrized form. 
To prepare the way for the study of the Klein-Gordon equation, we then briefly review 
some useful aspects of relativistic quantum mechanics in Section 3. 

Our main results are described in Sections 4 and 5 where we apply the formalism 
outlined above to the relativistic particle. In Section 4, we construct a position operator 
which commutes with the constraint. Its eigenstates are the Newton- Wigner states, and 
in fact the operator is essentially the same as the Newton- Wigner operator [22]. The 
associated probabilities on spacelike surfaces are those one would anticipate on the basis of 
the Schrodinger equation which is the square root of the Klein-Gordon equation. Another, 
different, candidate expression for the probability associated with a section of spacelike 
surface is the flux of the Klein-Gordon current (with the sign of the negative frequency 
part changed, as in Eq.(1.7)). This probability is related to a different set of position states 
which are non-orthogonal but relativistically invariant. There are, therefore, even at this 
simple level of canonical quantization, two distinct quantizations of the relativistic particle, 
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which are not equivalent (corresponding loosely speaking to "quantize then constrain" 
versus "constrain then quantize"). We shall refer to them as the Klein-Gordon (KG) and 
Newton- Wigner (NW) quantizations. 

In Section 5 we consider the decoherent histories analysis of the system in the KG 
quantization. We compute the decoherence functional for histories which cross a surface 
of constant x° in a spatial region A (or in its complement A). Since the paths move 
backwards and forwards in time, the notion of crossing is ambiguous and needs to be 
carefully defined. We show that the essentially unique notions of crossing associated with 
the KG quantization are first and last crossing. We thus obtain probabilities associated 
with the spacelike surface of the form of a Klein-Gordon inner product (with a sign change 
in the negative frequency sector). However, the histories are only approximately decoherent 
(with the off-diagonal terms proportional to the overlap of the relativistically invariant 
position states). 

The computation of the crossing probabilities hinges on resolving a subtle point: the 
path integral representation of the class operator for not crossing suggests that there is 
a non-zero amplitude that the particle will never cross a spacelike surface, contrary to 
intuition. This issue turns out in fact to be related to the problem of class operators 
which do not satisfy the constraint mentioned above. An important part of the analysis 
of Section 5 is a demonstration of how the class operators may be modified in a sensible 
way so that they do satisfy the constraint. The properly modified class operator for not 
crossing a spacelike surface then turns out to be zero, in agreement with physical intuition. 

In Section 6, we consider the decoherent histories analysis in the NW quantization. 
This is much simpler, being very similar in form to non-relativistic quantum mechanics. 
Decoherence is exact and the expect NW probability expressions are easily recovered. 

We summarize and conclude in Section 7. 
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This paper is related to a number of other works in the field. It exploits and extends 
the general formalism of the decoherent histories approach applied to quantum cosmology 
set out by Hartle [23], and more recently by Hartle and Marolf [16]. The application of 
this formalism appears to have been applied to particular models in only two other places. 
Whelan [24] considered probabilities on timelike surfaces for the relativistic particle (but 
not using the induced inner product structure, as here) . Craig and Hartle [25] have applied 
the formalism to a Bianchi IX quantum cosmological model. There is also some connection 
with the work on probabilities for non-trivial spacetime coarse grainings in non-relativistic 
quantum mechanics [26]. A more general investigation of these ideas applied to quantum 
cosmological models is currently being pursued in Ref.[27]. 



The very first testing ground for ideas about the quantization of reparametrization 
invariant systems is the parametrized non-relativistic particle. This is the usual non- 
relativistic particle but with the time coordinate t raised to the status of a dynamical 
variable, with conjugate momentum pt- Its action in Hamiltontian form is 



where a denote denotes differentiation with respect to the parameter s. N is a Lagrange 
multiplier enforcing the constraint 



where h is the usual Hamiltonian h = pj,/2m. Canonical quantization leads to the 
Schrodinger equation, 



2. THE PARAMETRIZED 



NON-RELATIVISTIC PARTICLE 




(2.1) 



H = p t + h = 



(2.2) 




(2.3) 



11 



In terms of dynamics nothing new is gained at this stage. But the interesting question is 
to see what the usual expressions for probabilities look like in the language introduced in 
Section 1. 

Following the general scheme, we normalize solutions to the constraint by first con- 
sidering eigenstates of H, as in Eq.(1.5). They are normalized using the auxiliary inner 
product according to 

*m)a = J dtdx V* Xk (x, t)9 x ,v(x, t) = 5(X - X')8(k - k') (2.4) 

Since H = pt + h, the solutions to the eigenvalue equation may be written 

*Afc(M) = — ^ e iXt i; k (x,t) (2.5) 
(2n)2 

where ip k (x,t) satisfies the Schrodinger equation. It follows that 



1 

2^ 



J dt J dx e- iXt+iX>t il>* k {x, t)4> k ,(x, t) = S(X - X')S(k - k') (2.6) 
The integral contains within it the usual inner product 

(V'jfe. Ms = J dx ipl(x, t)i/; k ,(x, t) (2.7) 

This has the important property that it is independent of time when the states obey the 
Schrodinger equation, so the time integral may be done in Eq.(2.6), pulling down a delta 
function S(X — X'), and it follows that 

^ k ,Ms = S(k-k') (2.8) 

This means that the expected Schrodinger inner product on surfaces of constant t follows 
from the induced inner product defined on the whole of spacetime. 

We may now ask for the probabilities in various situations of interest. Perhaps the 
simplest is the probability of finding the particle in the spatial region A at time tQ. In the 
usual approach it is of course 

p A = / dx |V(Mo)| 2 (2-9) 

Ja 
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To express this in the language of Section 1, we seek an operator which commutes with the 

constraint H and corresponds to the answer to "the value of x when t = to" • Following 
the general scheme this is 

X = (°° ds d ^fi- x(s) 6(t(s) - t Q ) (2.10) 
J-oo ds 

where x(s) and t(s) are the evolution of x and t using the constraint H as a Hamiltonian: 

x (s) = e lHs xe- lHs , t(s) = e lHs te- lHs (2.11) 

Since H = pt + h this is 

x(s) = e lhs xe~ lhs = x + ^ (2.12) 

m 

t( s ) = e tptS te- tptS = t + s (2.13) 

The integral over s may be done in Eq.(2.10) with the result 

X = x- P(t - t0) (2.14) 
m 

It is easy to confirm that this commutes with H. 

Since H and X commute, they possess a joint set of eigenstates u\%. The eigenvalue 
equation for X is 

Xu X x(x, t) = iit Ai (i, t) (2.15) 

with solutions 

u\dx, t) = —^e iXt g(x, t\x, t ) (2.16) 
(2n)2 

where g is the non-relativistic propagator. In the auxiliary inner product these are nor- 
malized according to 

(u\x, u X 'x')a = <5(A - A') 8(x - x) (2.17) 
The amplitude for an eigenstate of the constraint of the form 

V x ,(x,t) = -^e iXt i;(x,t) (2.18) 

(27T)2 
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to be in an eigenstate of X is 

{u\x^\>)a = t^t I dtdx e^ A - A V(M|MoMM) 
(2tt) J 

= S(X-X')i;(x,t ) (2.19) 

The probability is then computed from the expression 

/ dx (^ X "^\x)a (u x ^ x ,) A = 5(\" -X)5{\-\') [ dx\^t )\ 2 (2.20) 
J A J A 

Following the induced inner product prescription, we drop the delta-functions on the right, 

thereby obtaining the expect result for the probability, Eq.(2.9). 

The operator formalism with Eq.(1.8) allows one to ask a richer variety of questions 
than those normally considered in non-relativistic quantum mechanics. We may consider, 
for example, the question, "What is the value of t at a given value of x" ? The associated 
operator is 

^ mix — zn) 

T = t + ^ ^ (2.21) 

P 

A suitable operator ordering must be chosen, but the presence of the 1/p factor makes it 
difficult to turn this into a self-adjoint operator (see Ref.[28], for example). This operator 
arises in relation to the arrival time problem in non-relativistic quantum mechanics, an 
issue that has attracted a lot of recent attention in the literature [29]. 

Both of the above questions in non-relativistic quantum mechanics may also be anal- 
ysed using the decoherent histories approach. We will not go into the details here, except to 
make some simple observations that are related to the relativistic particle case we consider 
later. 

In the decoherent histories approach, the probability Eq.(2.9) may also be obtained 
using a standard non-relativistic path integral, in which one sums over paths which cross 
the surface t = to in the spatial range A. It is a property of this path integral that the 
paths cross this surface once and only once, and as a consequence of this, the histories are 
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exactly decoherent. As described in Section 1(C), exact decoherence of histories implies 
that records exist [7,21], or in other words, that the probability may be written in the form 
Tr(R a p) for some projection operator R a . This is thoroughly consistent with the existence 
of the self-adjoint operator (2.14) from which the probabilities Eq.(2.9) are derived in the 
operator approach. 

But now consider, by contrast, the probability for crossing a surface of constant x. 
In the decoherent histories analysis of this question (which is rather non-trivial [26]), 
the paths may cross the surface many times. Furthermore, it is found that the histories 
are typically not decoherent (unless an environment to produce decoherence is included, 
but we do not consider that case here), and this appears to be related to the multiple 
crossings. We cannot therefore deduce the existence of records and a probability of the 
form Tr(R a p). There would be an inconsistency with the operator approach here if there 
was a self-adjoint operator corresponding to this question. But interestingly, as we have 
seen, the corresponding operator Eq.(2.21) is not self-adjoint. The point therefore, is 
that multiple surface crossings and the associated lack of decoherence in the decoherent 
histories approach appear to be related to the absence of a self-adjoint operator in the 
operator approach. We will see more evidence of this in the case of the relativistic particle 
in Sections 5 and 7. 

3. GREEN FUNCTIONS OF THE 
KLEIN-GORDON EQUATION 

The Klein-Gordon equation has a variety of associated Green functions and it will be 
useful to briefly summarize them here. In order to agree with the notation of Ref.[30] 
(which we follow very closely), we use particle physics convention in which the signature 

of the metric is (H ). The positive and negative frequency Wightman functions G ± 
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are defined by 

G ±( x , y ) = * f dPk e -ik.(x-y) (31) 

1 ,y) (2tt) 3 J ko=±Wu 2u; k 1 J 

where cu^, = Vk 2 + m 2 . They satisfy the composition laws 

G ± = ±G ± oG ± , G ± oG T = (3.2) 

where here and in the remainder of this section o denotes the Klein-Gordon inner product 
(unless explicitly denoted otherwise). The causal Green function is defined by 

iG(x,y)=G+(x,y)-G-(x,y) (3.3) 

Its main property is that it propagates all solutions to the Klein-Gordon equation 

<j) = iGo<f) (3.4) 

It also obeys the composition law 

G = iGoG (3.5) 

The Hadamard function is defined by 

G^(x,y) = G+(x,y) + G-(x,y) (3.6) 

and obeys the composition laws, 

GW =iGo G^ = iG^ o G, G = o (3.7) 

All of the above are solutions to the Klein-Gordon equation. 
The Feynman Green function is 

iG F (x, y) = 6(x° - y°)G+(x, y) + 6(y° - x°)G-(x, y) (3.8) 

and satisfies 

(n + m 2 ^G F (x,y) = -S^\x-y) (3.9) 
16 



It obeys the composition laws 

Gp = iGp o Gp 



(3.10) 



Also of interest is the Newton- Wigner propagator 

G w (x, x°, y, y°) = * / rf 3 k e -<*-(*-») (3.11) 

which is the propagator associated with the positive square root of the Klein-Gordon 
equation 

= H (3.12) 

where h = \J — V 2 + m 2 . It is also useful to define a negative frequency Newton- Wigner 
propagator, given by (3.11) but with 1zq = —uj^, and this will be denoted Gjyw. It is easily 
seen that the Newton- Wigner propagator is related to the Wightman function by 

g nw (x, y) = 2i -Qj) G+ ( x i y) = - 2z ^o G+ ( x ' v) ( 3 - 13 ) 

Some of these Green functions can be obtained from a path integral of the form (1.15), 
(1.16). An unrestricted sum with T integrated over an infinite range yields the Hadamard 
function GW. (See Fig(3.1)). A half-infinite range, < T < oo, yields iGp, where 
Gp is the Feynman Green function. (See, for example, Ref.[31]). The Newton- Wigner 
propagator can also be obtained from (1.15), (1.16) by summing over all paths from y to 
x which never cross the surface of constant x°, except when they end at the point x. (See 
Fig. (3. 2)). More details of this construction are discussed in Sections 5 and 6. (See also 
Ref.[30]). 

From the path integral representations, one can see that G^ corresponds to the opera- 
tor 5(H), which is essentially the identity on the constraint subspace (and so we effectively 
have 5(H)\4>) = \<p) for solutions to the constraint). This is perhaps confusing since 
does not in fact propagate positive and negative frequency solutions to the Klein-Gordon 
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equation (it is the causal Green function G that does this job, via Eq.(3.4)). The reso- 
lution of this is the choice of inner product, does in fact propagate all solutions if 
they are attached with the induced inner product (1.7). For suppose we have a solution 
= 0+ + 0- . Then 

oj cj> = (G + + G~) oj (0+ + </>-) 
= G + o KG 0+ - G~ o KG 
= (G+-G-)°kg (<P + + <P~) 

= iGo KG <f> (3.14) 

In this sense, G^ is effectively equivalent to G. 

It is also interesting note in this connection that it was claimed in Ref. [30] that there is 
no path integral of the form (1.15), (1.16) that will yield the causal propagator G directly. 
Whilst this is still in some sense true, one can see that it depends on how the initial states 
are attached: the path integral for G^ but with initial states attached using the induced 
inner product does in fact effectively give the causal propagator G. 

We may now consider the form of the decoherence functional for the Klein-Gordon 
equation (we follow the construction of Ref. [16]). We take it to be of the form (1.14). We 
take a fixed pure initial state and sum over a complete set of final states. This gives 

D(a, a') = J^&f °I Ca oj V)% °I C a > oj </>)* (3.15) 

where note here that we use the induced inner product. (The normalization factor is unity 
in this case). Since ipj denotes a complete set of positive and negative frequency solutions, 
it is easy to show that 

Y,r f (x)iP f {y) = G^{x,y) (3.16) 
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Furthermore, since C a are solutions to the constraints, the action of changes nothing, 

GO-) 

/ — Ca? so we have 

D(a,a') = r°lCi,o I C a o I if; (3.17) 
Finally, when there is exact decoherence, D(a, a') = for a ^ a/, the probabilities are 

p(a) = D(a, a 1 ) = D(a, at) 

a 1 

= lP*o I C a o I iP (3.18) 

In the KG quantization, the induced inner product becomes the modified KG inner product 
(1.7), with initial states normalized in this inner product. In the NW quantization, states 
obey the Schrodinger equation, they are normalized in a Schrodinger inner product and o 
is taken to be that inner product in the decoherence fucntional. 



4. AN OPERATOR APPROACH FOR 
THE KLEIN-GORDON EQUATION 



We now describe the use of operator methods to obtained probabilities associated with 
the Klein-Gordon equation. The relativistic particle is described by the constraint 

H = pl-p 2 -m 2 = (4.1) 

where the canonical variables x^, p v obey the commutation relations 

\x», v v \ = -irfv (4.2) 

We are interested in the question, 'What is the value of x k when x° = r?". As indi- 
cated already, there are potentially many ways of formulating and answering this question. 
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We will first use the operator methods of Section 1(B). Following Eq.(1.8), the operator 
expressing this question is 

"oo i j„,0/ 



X k = J°°Js \{ d -^MAs)-r)}x\s) 
where { , } denotes the anticommutator, and we have 



(4.3) 



x°(s) = x +p°s (4.4) 
x\s) = x l +p l s (4.5) 

The object of interest is therefore given by 

X k = x k -^{\,(x Q -r)} (4.6) 
z p 

and it is easily seen that this commutes with H. One might anticipate that the factor 
may present problems in turning this into a self-adjoint operator, but this problem does 
not arise since we are looking for eigenstates of X % which also satisfy the constraint, and 
this bounds pq away from zero. (Essentially the same operator was also considered by 
Marolf [14]). 

We choose a momentum representation, in which 

,k ■ d J) ■ d 



and X k is 



x K -> -%—- , x v -> -i— J (4.7) 
dpk dp" 



Q pk q k k T 

X k = -l-z— + ^tttt - i nx0 + -7T 4.8 

1 dp k p°dp° 2(p ) 2 p° 1 



(recalling that p k = —p^ with our choice of signature). This is self-adjoint in the (momen- 
tum space version of) the auxiliary inner product Eq.(1.4). The eigenstates of X k are the 
functions 

f(p) = ^fj2^ h e^-P'* g(p ■ p) (4.9) 
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where g is any function of p -p, and the eigenvalue is x k . For these to be eigenstates of the 
constraint we also need to choose g = 8(p ■ p — m 2 ). Introducing the eigenstates \p) of p^, 
(where {p\p') = 5^\p — p')), the eigenstates of X 1 may be written 

|x, r) = J d 4 p (2p°)5 e iA-ip-x 6{p 2 _ m 2 } |p) 

(27T)3/2 7 ' P+) + (2^/2 i (2 ^i e IP ^ 

= |x,r+) + i|x,r-) (4.10) 

(The factor of z in the negative frequency term, not present in other definitions of these 
states [32,30], does not in fact make any difference.) Here, we have introduced the mo- 
mentum states |p±) on the positive and negative frequency sectors which are normalized 
according to 

(p±|p / ±) = 2^(p-p / ) (4.11) 

and also (p ± |p /= F) = 0. The states |x, r) are the Newton- Wigner states [22]. They are 
orthogonal at equal times, and satisfy the completeness relation 

1 = J d 3 x |x, r+) (x, t + I + J d 3 x |x, t-) (x, T - | (4.12) 

The probability of entering a region A at time r is given by 

p A = f d\ |(x,r+|^)| 2 + / d 3 x |(x,r-hA>| 2 (4.13) 
J A J A 

where the states ip±(~x., r) = (x, r ± \ip) obey the Klein-Gordon equation and its posi- 
tive/negative square root, and are normalized in a Schrodinger inner product. 

The Newton- Wigner states could also have been obtained by solving the constraint 
classically and then considering the eigenstates of the operators 

X k = x k ± - (4.14) 



a/p 2 + m 2 
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The Newton- Wigner states therefore correspond to "constraining before quantization" . It 
is also important to compare with the position operator introduced by Newton and Wigner 
[22], which is, on the surface x® = 0, in the momentum representation 

x ™ = -*ik-i v (4 ' 15) 

This is in fact the same as the operator 



(Fx |x, t+) x K (x, t + I (4.16) 

with r = 0, in terms of the Newton- Wigner states above. Eq.(4.15) is not the same as 
(4.6), since the constraint holds in Eq.(4.15), but has not yet been imposed in Eq.(4.6). 
Eq.(4.15) is, however, the same as Eq.(4.14), with r = 0, once one recognizes that the 
inner product structure (4.11) requires the replacement 

x k - - 2w p 2 i- r = -i—- - » 4.17 

d Pk (2w p )2 dPk 2^ 

There is therefore agreement with the earlier work of Newton and Wigner. 

An alternative way of defining position states is to first consider eigenstates of the 

position operator x^ on the auxiliary Hilbert space, and then project onto the constraint 

subspace using 5(H). This corresponds to quantizing before constraining, and yields 

I \ ^~ f dp zp'X I \ ^~ f dp zp'X I \ 

|X)= (2^4=. p 2V |P+)+ (2^i =_. p 2^ e |p - } 

= \x+) + \x-) (4.18) 

Unlike the Newton- Wigner states, these states are Lorenz-invariant. Furthermore, they 
not orthogonal, since 

(x\y) = G^(x,y) (4.19) 

although they are approximately orthogonal in the sense that G^(x,y) decays when x 
and y are separated by more than the Compton wavelength m _1 . They also obey a 
completeness relation 

l = i J d 3 x (j£+)cb(x + | - |x-)^b(x-|) (4.20) 
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These properties of the states \x) are reminiscent of the coherent states, and suggest that 
the probability for crossing a spacelike surface x® = r in the region A may be taken to be 



The states <j>±(x) = (x± |0) are positive/negative frequency solutions to the Klein-Gordon 
equation. The minus sign in the second term ensures that the expression is positive, and 
in the limit A = U-fi this expression becomes the norm of <p m the induced inner product, 
as required. 



We now come to the main point of this paper which is to use the decoherent histories 
approach to compute the answer to the question, "What is the probability that the particle 
is found in the spatial region A at time x° = r?". In this Section, we consider the KG 
quantization, with the aim of obtaining Eq.(4.21), and we consider the NW quantization 
and Eq.(4.13) in Section 6. 

The decoherence functional is given by Eq.(3.17). We take a pure initial state, and we 
sum over a complete set of positive and negative solutions in the final state. The main 
issue is to compute the class operators C a (x" \x') corresponding to crossing x° = r either 
inside the region A or outside it, in its complement A. We expect that these operators 
can be obtained by a sum over paths which either cross or do not cross the region (in a 
sense to be made more precise below). However, in order to tackle a subtlety first noted 
in Section 1 concerning the definition of the class operators, we first need to consider a 
simpler question. 




(4.21) 



5. DECOHERENT HISTORIES ANALYSIS 



IN THE KLEIN-GORDON QUANTIZATION 
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5(A). The Class Operator for 



Not Crossing a Spacelike Surface 



In computing a class operator of the form C a (x", x') we sum over paths from x' to x" 
satisfying some condition specified by the coarse graining a. However, this construction 
appears to allow for the possibility of defining a coarse graining consisting of paths from 
x' to x" which never cross a given spacelike surface. (See Fig. (5.1)). This is at first sight 
disconcerting. Every classical trajectory of the relativistic particle crosses every spacelike 
surface (unless it is tachyonic). And in the quantum theory, a solution to the Klein-Gordon 
equation cannot be zero on one side of a spacelike surface and non-zero on the other. 

To investigate this issue, we consider the following coarse graining: paths from x 1 to 
x" which either cross or do not cross the surface x° = r. That is, we are interested in 
the question, given a solution to the Klein-Gordon equation, what is the probability that 
the particle will be found on the spacelike surface x® = r, or will never be found on that 
surface? Clearly the answer to this question must be probability unity for crossing the 
surface, and probability zero for not crossing. 

The sum over paths which do not cross x° = r is easily constructed using the path 
integral representations^. 15), (1.16). We take T to have an infinite range and we sum 
over all paths which never cross the surface. In Ref.[30] it was shown that this is in fact 
equivalent to a method of images construction for g(x" , T\x', 0), and we obtain, for the 
class operator for paths restricted to not cross the surface, 





(5.1) 



where x denotes the reflection of the point x about x = r, that is, 



x» = (2r-x u ,x) 



(5.2) 
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C r vanishes on x ' = r and x°" = r. Both G^(x",x') and G^(x",x') are solutions to 
the KG equation, but the presence of the 6>-functions means that C r does not satisfy it, 

(u + m 2 )C r (x",x') = 25(x°" -t) e(x ' - r) d G^\T,x"\x \x') (5.3) 

where e(x) is the signum function. This is the difficulty with the path integral-defined 
class operators mentioned in Section 1. The sum over paths which always cross x° = r, 
which we denote C c {x" \x'), is constructed using the relation 

G {l \x", x') = C c (x", x') + C r (x", x') (5.4) 

from which one can see that C c (x", x') will also not satisfy the constraint (although to 
actually compute the crossing class operator we will use a different method below). 

Since the constraint is associated with reparametrization invariance, it might seem that 
failure to satisfy it is associated with a breaking of that invariance. However, the connection 
between constraints and invariance can be rather subtle for the case of reparametrizations. 
In particular, the path integral (1.15), (1.16) between fixed end-points over paths restricted 
to pass through certain spacetime regions is clearly reparametrization invariant. These 
issues are discussed more fully in Refs.[l,33]. 

Note also that we have now encountered two potential difficulties which are not obvi- 
ously related: firstly, the possibility of a physically unreasonable result, and secondly, the 
fact that the naturally defined class operator does not satisfy the constraint. We will now 
see that both of these problems are solved simultaneously by appropriate modification of 
the class operator. 

Following the suggestion of Hartle and Marolf [16], we will deal with this issue by 
replacing C r by another object C' r which satisfies the most important boundary conditions 
defining C r but which also satisfies the constraint everywhere. The boundary conditions 
satisfied by C r are that, (a), it vanishes on x®" = r, x®' = r and when x" and x' are on 
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opposite sides of the surface. One might also be tempted to impose that, (b), C' r coincides 
with (the non-zero) C r when x" and x' are on the same side of the surface. It appears 
to be impossible for a function satisfying the constraint everywhere to satisfy all of these 
conditions. The essence of the non-crossing propagator appears to be contained in the 
conditions (a), so we drop conditions (b). The unique solution to the constraint equation 
satisfying conditions (a) is then quite simply 

C' r (x",x') = (5.5) 

This is because for fixed y, C' r (x", x') is zero for all values of x" on the opposite side of the 
surface to x' . This means that both C' r and its normal derivative are zero on all spacelike 
surfaces on the opposite side of the surface x® = r, and the solution to the Klein-Gordon 
equation with these conditions is simply C r = 0. The only way of getting a non-zero 
results as x" moves from the opposite side to the same side as x' is to have a discontinuity 
in the normal derivative, but this can only be achieved with a delta-function source, as in 
Eq.(5.1). 

The conclusion (5.5) for the modified class operator for non-crossing paths implies that 
the modified class operator for paths that always cross is quite simply G^ l \x,y). It now 
follows that histories partitioned according to whether or not they cross x Q = r are exactly 
decoherent, the probability of not crossing is zero, and the probability of crossing is 1. We 
have therefore shown that by sufficiently careful treatment of the class operators and their 
boundary conditions, we obtain the expected and physically sensible result. 

5(B). Crossing Propagators and 
the Path Decomposition Expansion 

Having resolved the issue of how to modify class operators that do not satisfy the 
constraint, we may now turn to the issue of computing the class operators for crossing the 
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spacelike surface i - r in a spatial region A. These operators (before modification) will 
be constructed by a path integral of the form (1.15), (1.16) in which the paths cross the 
spacelike surface in the region A. Because the paths go backwards and forwards in time, 
they will typically cross a given spacelike surface many times, so the notion of crossing 
needs to be specified more precisely. We will see however, that for the path integral 
representations of the Klein-Gordon propagators, first and last crossings are in fact the 
only useful notions of crossing. 

A very useful result for our purposes is the path decomposition expansion, or PDX 
[34,30]. For a propagator of the non-relativistic form (1.16), it implies that when x" and 
x' are on opposite sides of the surface of constant a; , we have 

g(x", T\x', 0) = 2% e(r — x° ) J dt c J (fix g(x", T\x, t c ) dog{x, t c \x' , 0) (5.6) 

This formula is obtained by partitioning the paths in the sum over histories (1.16) according 
to the parameter time t c and position x at which they cross the surface of constant x® for 
the first timet. (See Fig. (5. 2)). If we partition according to the last crossing, we get 

g{x", T\x', 0) = -2i e(x°" - r) J dt c J d 3 x g(x", T\x, t c )dog(x, t c \x' , 0) (5.7) 

When x" and x' are on the same side of the surface, there is the possibility of paths between 
these points which do not cross the surface. The appropriate formulae then are 

g(x",T\x',0) = g r (x",T\x',0) 

+ 2ie(r-x 0/ ) J dt c J d 3 x g(x", T\x, t c ) dog{x, t c \x', 0) (5.8) 

for the first crossing and 

g(x",T\x'^)= g r (x",T\x',0) 

-2ie(x°"-r) J dt c J d 3 x g(x", T\x, t c ) dog{x, t c \x', 0) (5.9) 

t The full version of the PDX actually involves the normal derivative of the restricted prop- 
agator g r , but in this simple case, g r may be computed using the method of images and it 
follows that dog r = 2d$g, which is what is used in Eq.(5.6). See Ref.[30] for more details. 
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for the last crossing, where g r denotes the restricted propagator given by a sum over paths 
which never cross the surface. The formulae (5.6)-(5.9) imply that the sum over paths 
which cross the surface is given by either (5.6) or (5.7), irrespective of whether x' and x" 
are on the same side or opposite sides of the surface. (But only in the latter case are these 
expressions then equal to the full propagator). 

These results were used in Ref. [30] to derive the composition laws of relativistic propa- 
gators from the path integral. Here, we note that the propagators for first or last crossing 
the surface x® = r in the region A are readily obtained by simply restricting the d^x 
integration to the region A. 

5(C). First and Last Crossing Relativistic Propagators 

Turning now to the relativistic propagators, the class operator for crossing the surface 
x® = t in the region A is 

/oo 
dTg A (x",T\a/,0) (5.10) 
-oo 

We first take the case where g&{x" ', T, x', 0) is the sum over paths from x' to x" in fixed 
proper time T which cross the surface for the first time in A, that is, 

g A (x",T\x',0) = 2i e(r - x ') / dt c d 3 x g(x", T\x, t c ) d g(x, t c \x' , 0) (5.11) 

JO J A 

(See Fig. (5. 3)). This is valid for x" , x' on either the same or opposite sides of the spacelike 
surface. Inserting in Eq.(5.10), writing the integral as a sum of two parts corresponding 
to the positive and negative ranges of T, and changing variables to v = T — t c , u = t (see 
Ref. [30] for more details), this yields, 

C f A (x", x') = -2i e(r - x°) [ d 3 x 

J A 

(5.12) 
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G F {x",x)d G F (x,x') - G F {x",x)d G F {x,x') 



This is the formula for first crossing the region A for all end-points. It is also conveniently 
written, 

C{(x",x') = - [ d 3 x \G^(x",x)~doG(x,x') + e(x°" - r)e(r - x ') G(x" , x)doG^\x, x') 
J A L 

(5.13) 

It is also of interest to consider the class operator defined by the last crossing, which is 
easily shown to be 



C e A (x",x') =2i e(x v " -t) / d 6 x G F (x",x)d G F (x,x') - G F {x" ,x)d Q G F {x,x') 



l\ t '" 



1* („ J\ 



L 



= / d 6 x 



e(x°" - r)e(r - x°')G^\x", x]d Q G{x, x') 
+ G(x",x)doG^\x,x') 



(5.14) 



When the initial and final points are on opposite sides of the surface, we have 



e(x o" _ r)e(r _ x o' } = 1 



(5.15) 



It is then convenient to average the first and last crossing class operators to obtain 

C A (x", x') = i (C{(x", x') + Ci(x", x')) 

= -- [ d 3 x \g( 1 \x",x)W G(x,x')+ G(x",x)WoG^\x,x') 
2 J A 1 J 

= ij d 3 x {G + {x\x)%G + {x,x')-G-{x'',x)%G-{x,x')^ (5.16) 

It is then readily confirmed, using the properties Eq.(3.2) and (3.7), that this class operator 
become in the limit that A becomes IB? , as expected. 

As one of x" or x' is moved from the opposite to the same side of the surface, the class 
operator undergoes a discontinuity. This is reflected in the fact that it does not satisfy the 
constraint. The first crossing class operator, for example, satisfies the equation, 

(a" + m 2 ) C{(x", x') = 2 e(r - x ') 5(x°" - T)d G^(r, x"\x ' , x') (5.17) 
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when x" is in A and zero otherwise. Note that Eq.(5.3) and Eq.(5.17) are consistent, since 
C/± + C r = when A = iR 3 , and GO) satisfies the constraint. As in Section 5(A), 
some doctoring of this basic class operator must therefore be carried out before we get 
the final expression, for a modified class operator C^, which satisfies the constraint. It is, 
however, clear in this case how to proceed. From the above that the obvious candidate 
is to take C' A to be given by Eq.(5.16) for all values of the end-points x" , x' , whether 
they lie on the same side of the surface or opposite sides. This is clearly a solution to the 
constraint everywhere. It matches the path integral-defined object when the end-points 
are on opposite sides of the surface. Furthermore, when A = iR 3 , it is equal to G^ 1 ) , so is 
consistent with the modified class operator for not crossing. We wil return below to the 
question of a more general prescription for constructing modified class operators. 

We may now consider the decoherence functional for histories which cross the spacelike 
surface x® = r either in the region A or in its complement A, for an initial state ip = 
ip+ +ip~. The off-diagonal terms of the decoherence functional are 

D(A, A) = r oj (tf A )t 0/ C' A oj i, (5.18) 

We have 

C' A oj^ = i J d 3 x (G + {x",x)^ + {x) + G-{x",x)%i)-{x)^ (5.19) 
so the decoherence functional is 

D(A, A) = J d 3 x Sy (r+(y)%G+(y, x)t^+{x) + V* (y)^G _ (y, 

(5.20) 

The key feature of this expression is that the decoherence functional is not exactly diagonal. 
It is, however, approximately diagonal in the sense that the two-point functions G^(y,x) 
decay for increasing spatial separations. In particular, we expect that approximate di- 
agonality can be obtained if both regions A and A are much larger than the Compton 
wavelength m" 1 (the decay length scale of G ± {x : y)). 
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Given decoherence, the probability for crossing A is then 

V>* °I C' A oj xj, = i J d 3 x (^* + (x)%^+{x) (5-21) 

This is exactly the expected answer, coinciding with Eq.(4.21), although recall that the 
probabilities are only approximately defined because of approximate decoherence. 

5(D). A General Prescription for 
Constructing Modified Class Operators 

We have so far constructed the modified class operators using some general arguments, 
but the question remains as to whether it is possible to find a more general formula for 
constructing them. Connected to this is the question of how the modified class operators 
are related to the original path operators, such as Eq.(5.10), which were defined using path 
integrals in a simple and obvious way. 

Consider first, therefore, the question of why the expressions (5.10) and (5.11) fail to 
satisfy the constraint. Using the fact that the propagators of the form g(x, t\x , 0) satisfy 
the Schrodinger equation, it is easy to see that (5.10) fails to satisfy the constraint at 
x" because of the finite integration range for t c . Recall that t c is the parameter time of 
first crossing and is not a physically observable quantity. Because it is unobservable, and 
because the total parameter time T is integrated over an infinite range, it seems reasonable 
to explore the possibility that t c could also be integrated over an infinite range, in such a 
way that a solution to the constraint equation is obtained. 

Proceeding along these lines, one can see that one way to obtain a solution to the 
constraint is to extend the integration range of t c to — oo < t c < oo, but with a signum 
function e(t c ) included. That is, we define the modified class operator 

f roo poo r — i, 

C' A (x", x') = 2i e(r - x° ) dT dt c e(t c ) / d 3 x g(x", T\x, t c ) d g(x, t c \x' , 0) 

J —oo J —oo J A 

(5.22) 
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Now we note that 

/oo 
dt c e(<c) g(x, t c \x', 0) = iGp(x, x') + iG* F (x, x') 
-oo 

= ie(r - x ') G(x, x') (5.23) 

(where recall x = r). Furthermore, we have that 

9q (e(r - x ') G(x, x')) = e(r - x ') ^,G(x, x') (5.24) 

since G(x, x') vanishes when r = x '. The modified class operator is therefore 

C' A {x",x') = -2 f d 3 x G^(x",x)^ G(x,x') (5.25) 
J A 

which, apart from the factor of 2, coincides with the first term in Eq.(5.13). It is easy to 
see that the second term in Eq.(5.13), but without the e factors (and again up to a factor 
of 2) may be obtained by a slightly different modification of the integration range over 
parameter times, that is, 

f roc roo r — y 

2ie(r-x°) / dT e(T - t c ) \ dt c / d 3 x g(x", T\x, t c ) d g(x, t c \x, 0) 
J— oo J —oo J A 

= -2 [ d 3 x G{x",x)dQG^\x,x') (5.26) 
J A 

Hence averaging these two results gives the class operator Eq.(5.13), but crucially, without 
the e factors that cause (5.13) to fail to satisfy the constraint. We may finally perform a 
further averaging with the last crossing versions of the above modified class operators to 
obtained a modified class operator with the following properties: it satisfies the constraint 
everywhere, and, when x" and x' are on opposite sides of the surface, coincides with the 
path-integral defined object (5.16). 

In summary, we therefore now have two different methods of defining modified class 
operators which satisfy the constraint. One is to use the usual path integral construction 
to compute the class operators when x" and x' are on opposite sides of the surface, and 
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then declare that this is valid for all values of x" , x' . The second is to use the path integral 
defined object but modify the integrations over the unphysical parameter time labels, as in 
Eqs.(5.22), (5.26). The two methods are equivalent for the simple model of this paper. The 
benefit of introducing the second method is that it shows that the modification procedure 
does not fundamentally modify the class of paths in configuration space summed over in 
the path integral, only the way their parametrizations are summed over. The modified 
class operators therefore still have equal claim to be a sum over paths which cross the 
surface. 

Theses issues concerning modified class operators will be taken up in more detail in 



The results of Section 5(C) show that it is a partition of paths according to their first 
and last crossing of a spacelike surface that leads to the expected Klein-Gordon probability 
expression Eq.(4.21). It is, however, of interest to explore other notions of surface crossings. 
This is partly by way of a digression, but it is also relevant to the recovery of the Newton- 
Wigner probability expression. 

We begin by considering the first-crossing path decomposition expansion Eq.(5.6), 
with T integrated over a half-infinite range. We restrict to the case x Q " > r > x ', and for 
simplicity restrict to the positive frequency sector. In this subsection we will take A = iR 3 . 
We thus obtain 



(averaged with the last crossing PDX this gives the composition law in Eq.(3.2)). From 
Section 3 we also have that 



Ref.[27]. 



5(E). A Multiple Crossings Decomposition? 




(5.27) 



2id Q G + ( 



x", x) = Gnw(x", x) 



(5.28) 
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and Eq.(5.27) becomes 

G + {x",x') = J d d x G + {x",x)G NW {x,x') (5.29) 

In this expression the Newton- Wigner propagator therefore represents paths that start at 
x , move forwards and backwards in time but without crossing the final surface, except to 
end on it at point x. (See Fig. (3.2)). Similarly, one can see from the last crossing PDX 
Eq.(5.7), that a sum over paths which move backwards and forwards in time but without 
crossing the initial surface, except to start on it, is — 2iG + (x", x') 8q. (See Fig. (5. 4)). This 
is in fact again the Newton- Wigner propagator Gjyw(x, x'), as may be seen from Eq.(3.13), 
and we may write 

G + {x",x') = J d z x G NW (x",x)G + (x,x') (5.30) 
Using these two relations, we may carry out an iteration of Eq.(5.29) to yield 

G+(x", x') = J d 3 Xl d 3 x 2 G NW (x", x 2 )G+(x 2 , xi)G NW { Xl , x') (5.31) 

The two Newton- Wigner propagators represent restricted propagation to the first crossing 
and from the last crossing of the surface, and the propagator G + {x2- l x\) is unrestricted 
propagation between two points which both lie on the spacelike surface. (See Fig. (5.5)). 

It is now reasonable to ask whether this expression can be further decomposed accord- 
ing to the detailed number of surface crossings entailed in the path integral representation of 
G + (x2, x\). Interestingly, this does not appear to be possible. For suppose we apply a first 
or last crossing expansion of the type Eq.(5.29) or Eq.(5.30) to the propagator G + (x2, x\). 
The point is that all three points involved (initial point, final point, crossing point) all have 
the same value of x®, hence the expression would contain the Newton- Wigner propagator 
Gnw{ x iV) at x ° = 2/°- But this is simply the delta function #( 3 )(x — y), leading to a 
trivial result. Hence a further decomposition according to the specific number of crossings 
appears to be impossible. 
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The explanation for this is that in a path integral representation of G + (x2, x\), generic 
paths cross the surface an infinite number of times, and the set of paths crossing a finite 
number of times is of measure zero. This result is due in essence to Hartle [23], who 
considered a lattice version of the Euclideanized sum over histories. He attempted to factor 
the usual propagator of non-relativistic quantum mechanics across an arbitrary surface in 
spacetime by partitioning according to the number of crossings each path makes. He 
showed that paths with a finite number of crossings generally have zero amplitude in the 
continuum limit, and deduced that such a factoring is not in fact possible. (This is not 
at variance with the path decomposition expansion, Eq.(5.6), which partitions the paths 
according to their first crossing). 

On the face of it, therefore, it might seem like there are a number of different notions 
of surface crossing. The above results show, however, that first and last crossings are the 
only ones that can be defined in this case, hence are the only useful ones for defining the 
class operators of interest here. 

6. THE NEWTON- WIGNER CASE 

Consider now the question of how to obtain the Newton- Wigner probability Eq.(4.13) 
from the decoherent histories approach. Since this is like non-relativistic quantum me- 
chanics but with the Hamiltonian h = \J — V 2 + m 2 , it is simpler than the previous case 
and we describe it only briefly. 

The decoherence functional is given by Eq.(3.17) in which the inner product o is the 
Schrodinger inner product, with the initial states normalized in this product. The main 
issue is the construction of the class operator representing crossing a spacelike surface 
x° = r in a spatial region A. We first give the result and then explain its origin. (For 
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simplicity we concentrate on the positive frequency sector only). It is clear that the class 
operator is 

Ca(x",x')= / (fix GjvviK^", x, t) Gnw(x.,t,x') (6.1) 

Ja 

Note that when A = IR? this gives the standard composition property of the Newton- 
Wigner propagator 

Gnw(x",x) = J (fix Gnw{x" ,?l,t)Gnw{il,t,x) (6.2) 

Inserting in the decoherence functional (3.17), it is readily shown this gives exact deco- 
herence, and the probabilities coincide with Eq.(4.13) (this is very similar to standard 
calculations in non-relativistic quantum mechanics). It is necessary also to use here the 
fact that the Newton- Wigner propagator is the overlap of two Newton- Wigner states, 

G NW (x,x \y,y ) = (x,x°\y,y } (6.3) 

There are (at least) two path integral representations of the path integral for the 
Newton- Wigner propagator that lead to the class operator (6.1). The first is the one of 
the standard non-relativistic form: 



G(x", t", x', t') = j VxVp exp ^ jT dx° p • ^ - yW™ 2 



(6.4) 



(The configuration space form of this path integral may also be considered, but the measure 
is then rather complicated [32]). In this representation, the paths move forwards in the 
time coordinate x°. Summing over paths from x' to x" which pass through A on an 
intermediate spacelike surface then yields the class operator (6.1). 

A second and perhaps more interesting path integral representation of Gj^\y( x '' '> x ') ^ s 
the one mentioned in the discussion of the path decomposition expansion of the previous 
section. This is to use a path integral representation of the form (1.15), (1.16), in which 
the paths summed over do not cross the final surface except to end on it at x", depicted in 
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Fig. (3.2). (See also Refs. [23,30]). Equivalently, they can be restricted so that they start at 
the initial point x' but do not cross it thereafter (see Fig. (5.4)). In fact, it is easy to show 
from these representations, using Eq.(6.2), that Gnw(x", x') is obtained more generally 
by choosing any surface of constant x° lying between initial and final points, and then 
summing over paths which cross it once and only once, as depicted in Fig. (5.6). The first 
two representations then correspond to the limit in which the intermediate surface tends 
to the initial or final surface. From this third representation, we see that the class operator 
Eq.(6.1) is obtained by summing over paths which cross the intermediate spacelike surface 
once and only once, in the spatial region A. 

Note that there is no conflict here with the statement in Section 5(D) that the set 
of paths crossing a surface only a finite number of times is of zero measure in the set 
of all paths. Section 5(D) concerned path integral representations of the Klein-Gordon 
propagators, which involve a sum over all paths between two points, and the set of paths 
making single crossings of given surface are indeed insignificant. However, the path integral 
representation of the Newton- Wigner propagator considered is defined from the outset by 
a sum over the much smaller class of paths which cross an intermediate surface only once. 



7. SUMMARY AND DISCUSSION 



This paper is a first step in a programme whose general aim is to supply a reasonable 
predictive framework for quantum cosmological models. In connection with that aim, 
the main achievement of this paper is the derivation of the probability formula Eq.(5.21) 
(or Eq.(4.21)) from the decoherent histories approach, and the demonstration that the 
associated histories are approximately decoherent. 

Along the way we derived a number of other relevant results. We showed how to 

37 



modify in a physically sensible way class operators which do not satisfy the constraints. 
We also showed that first and last crossings are essentially the only ways of defining surface 
crossings (in the Klein-Gordon quantization). In particular, partitions of paths according 
to multiple crossings are not possible. These results will be relevant to more complicated 
models in quantum cosmology. 

The main result of Section 4 is the computation of an evolving constraints operator for 
the relativistic particle, proof that its eigenstates are the Newton- Wigner states, and that 
the operator is essentially the same as the Newton- Wigner operator. We discussed some 
novel path integral representations of the NW propagator in Section 6, involving single 
surface crossings and computed the decoherence functional (although noted that it is very 
similar to the case of non- relativistic quantum mechanics). 

We do not expect that a Newton- Wigner quantization based on a Schrodinger equation 
such as Eq.(3.12) will be relevant to more complicated models in quantum cosmology, since 
it is only under very special circumstances that the constraint may be solved to produce 
a real, positive Hamiltonian h to go in Eq.(3.12). The comparison with this case has, 
however, proved quite useful in the present paper. Furthermore, here, our starting point 
for an operator quantization was the evolving constants method, based on Eq.(1.8), which 
does not require the solution to the constraints, and this method will be valid for more 
complicated models (and indeed has been used already in such a context [14,15]). We also 
note that we find agreement between the operator methods of Section 4 and the decoherent 
histories results of Sections 5 and 6. 

We may now return to the discussion initiated at the end of Section 2, on the rela- 
tionship between decoherence, surface crossings and the existence of self-adjoint operators. 
We see further evidence for this. It is striking that the NW quantization, which involves 
single surface crossings in the decoherent histories approach, yields exact decoherence of 
histories, whereas the KG quantization, which has multiple surface crossings, exhibits only 
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approximate decoherence. It is reasonable to conclude from this that the approximate 
nature of the decoherence is related to the fact that the paths in KG quantization go back- 
wards and fowards in time and cross a surface many times. In fact, generally speaking, one 
might expect such multiple crossings to destroy decoherence altogether, since the paths 
may pass through both A and its complement A, but with single crossings they may pass 
through only one or the other. The interesting question is therefore why even approximate 
decoherence is obtained. For the free relativistic particle considered here, the answer is 
that the dominant contribution to the path integral representation of the class operators 
comes from the immediate neighbourhood of the classical path from x' to x", and this 
crosses the surface only once. Paths with multiple crossings therefore presumably belong 
to the quantum fluctuations about the classical path, and these may be neglected at suf- 
ficiently coarse-grained scales. Note also that the anticipated connection with self-adjoint 
operators holds up. The NW probability is associated with a self-adjoint operator, whilst 
the KG probability is not. 

These issues concerning multiple crossings become more complicated in non-trivial 
quantum cosmological models where typically even the classical paths have multiple surface 
crossings. This will be pursued elsewhere [27]. 
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Figure 3.1 The propagators ) and Gp(x",x') are obtained by a path integral 

of the form (1.15), (1.16), involving a sum over all paths from x to x (hence the paths 
may move both backwards and forwards in the time coordinate x°). An infinite range for 
i gives , x') and a half-infinite positive range gives iGp(x", x'). 
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Figure 3.2 The Newton- Wigner propagator Gnw{x" , x') may be obtained by a path 
integral of the form (1.15), (1.16), in which the paths may move backwards and forwards 
in time, with the restriction that they do not cross the final surface, except to end on it 
at the final point x" . 
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x' 

Figure 5.1 The class operator C r (x",x') for not crossing the spacelike surface x° = r is 
given by a sum over paths from x to x which never cross the surface. It is equivalent to 
a method of images construction. 




x' 

Figure 5.2 The path decomposition expansion (PDX) for the surface = r. A sum over 

paths from points x' to x" on opposite sides of a surface may be partitioned according to 

the position x and the parameter time at which it makes its first crossing of the surface. 
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x' 

Figure 5.3 The class operator C/^(x", x') for a first crossing of the surface x° = r in the 
spatial region A is obtained by summing over paths which cross the surface for the first 
time in the region A. 




Figure 5.4 A second representation of the Newton- Wigner propagator G]y\y(x" ', x') may 
be obtained by a path integral of the form (1.15), (1.16), in which the paths may move 
backwards and forwards in time, with the restriction that they do not cross the initial 
surface, except to start on it at the initial point x' . 
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Figure 5.5 The path integral representation of Eq.(5.31). The paths go from the initial 
point x' to their first crossing of the surface at x\, and this is represented by Gj\mr(x\, x'). 
The propagation from x\ to x 2 is unrestricted and is represented by G + (x\, X2). The paths 
make their last crossing at X2 moving to their final point x" , and this is represented by 
G NW (x",x 2 ). 




Figure 5.6 A more general path integral representation of the Newton- Wigner propagator 
Gnw(x" ,x'). In the path integral expressions (1.15), (1.16), the paths move backwards 
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and forwards in time with the restriction that they may cross a prescribed intermediate 
surface once and only once. The representation is independent of the choice of intermediate 
surface. The two previous representations (portrayed in Fig. (3.2) and Fig. (5.4)) in the limit 
that the intermediate surface tends to the final or initial surface. 
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